Comment on "Spin-1 aggregation model in one dimension" 
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M. Girardi and W. Figueiredo have proposed a simple model of aggregation in one dimension to 
mimic the self-assembly of amphiphiles in aqueous solution [Phys. Rev. E 62, 8344 (2000)]. We 
point out that interesting results can be obtained if a different set of interactions is considered, 
instead of their choice (the s = 1 Ising model). 

PACS numbers: 82.60.Lf,64.75.+g,64.00.Ht 



M. Girardi and W. Figueiredo W (henceforth GF) have 
considered the s = 1 Ising model in one dimension as 
a system that could possibly show some of the features 
of micellar systems, specially the existence of a critical 
micelle concentration (CMC). Several definitions of the 
CMC can be provided, since this phenomenon is actually 
a crossover, not a well defined phase transition. Follow- 
ing the classic work of H. Wennerstrom and B. Lindman 
|^| , they focus on the distribution of aggregates with n 
molecules, p(n). In real systems the curve np(n) is ex- 
pected to show a sharp maximum at n — 1, an interme- 
diate minimum and another maximum at some n corre- 
sponding to the micelles, typically of the order n « 10 2 . 
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where the arrows represent the +1 and —1 spins and 
the circle the one. Consider, instead, the following 
choice, in which we choose a left-right picture of the spins, 
instead of the usual up-down: 
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This is indeed the case, at least in some regimes, in the 
two- and three-dimensional models previously studied by 
the same group ||, f|, of which GF is a logical extension. 
On those works the difference between the minimum and 
the maximum at n was used as a sort of "order parame- 
ter" , whose vanishing could be associated with the CMC. 
In the case studied in GF, it was found that this func- 
tion is never bimodal and so a CMC cannot be properly 
defined. 

While their methods are not to be questioned, since 
they have carried out Monte Carlo simulations and, most 
importantly, it is possible to obtain an exact solution to 
the model, we feel that the model considered is not the 
best one if one wants to compare with real amphiphilic 
systems and that the failure to find a clear CMC seems 
natural in retrospect. As discussed in |^], the mean fea- 
ture of a micellar aggregate is the existence of a preferred 
aggregate size, which in real systems is a consequence of 
the effective geometry of the molecules (as exemplified 
by the well know picture of Israelachvili's |6|). A model, 
even a highly simplified one, should include some kind of 
preferred aggregate size in order to present a CMC. 

In a nutshell, the interaction considered in GF can be 
described by the following choice of spin interactions: 



We can expect this choice (which is basically an applica- 
tion of the potential used in Q to a lattice) to provide 

"micelles" consisting of spin pairs with orientation < > 

on a "solvent" of o spins, since this is the configuration 
energetically favored. This is indeed the case: e.g., we 
have checked that the slope of the chemical potential as 
a function of the logarithm of the amphiphilic density 
changes from 1 to a value 1/2, and this crossover corre- 
sponds to the CMC H ; this change in slope is not found 
with the GF model. 

However, this model is too simple to address the be- 
havior of np(n), since there is no possible aggregate be- 
tween n = 1 and n — 2 where the distribution would have 
a minimum. Our model can then be trivially extended 
to another one with to spins (to — 1 orientations for the 
amphiphiles and one for the solvent) Q]: 
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In Figure 1 we show the chemical potential versus the 
logarithm of the amphiphile density for the particular 
choice to = 5 to show the expected change in slope for 
1 to 1/4. We also plot include the corresponding values 
of np(n). It is clearly seen how the feature sought in 
GF is indeed obtained at ln(/3) « —18: at this point a 
unimodal distribution with a maximum at n = 1 turns 
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FIG. 1: Concentration of sites belonging to aggregates with 
different n (given in the legend) as a function of amphiphilic 



ference between the n = 1 and the n maxima) should 
perhaps be a better choice that the one in GF and pre- 
vious works (Refs. || (§)• 

We can also fix the concentration at some value and 
vary the interaction parameter (3 J. In Figure 2 we show 
results for p — 0.01. The results are quite similar to 
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FIG. 2: Same as Fig. 1 but now for a fixed amphiphilic density 
p = 0.01 and varying [3 J . 

into a bimodal with a second maximum at n — 4 (in 
general, one finds n — to— 1); the order parameter defined 
in GF is seen to vanish linearly. 

However, this point is quite far from the CMC one 
would define from the chemical potential, ln(p) « —6; on 
the other hand, the crossing between the n = 1 and the 
n = 4 lines does lie in this range. We suggest that this 
kind of criterion for the CMC (in general, np(n) = p(l), 
so that an order parameter could be defined as the dif- 



FIG. 3: Location of the CMC in the ln(p)-{3J plane, defined 
as p(l) = np(n) (solid lines) and as np(n) = (n — l)p(n — 1) 
(dashed lines), both for m — 5 and m — 10. 



the ones in Fig. 1, except that the chemical potential 
is now determined by our choice of density and does 
not show any particular feature. These two choices of 
crossing the CMC line yield the same results, which we 
have combined in Figure 3. There we present results 
for ttl — 5 and also for to = 10. We can see the CMC 
lines with the hp(n) = p{\) criterion show the expected 
dependence on m, with the slopes approaching a slope 
to — 2 at high interactions and low densities (since each 
aggregate contains to — 2 favorable bonds). The other cri- 
terion, based on the bimodality of the distribution (i.e., 
np(n) = (n — l)p(n — 1)) is seen to provide very differ- 
ent results, with a limiting slope of 1, independent on 
to. Nevertheless, the discrepancies between these two 
criteria can be expected to be smaller in real amphiphilic 
systems, since the energy will quickly grow for aggregates 
either smaller or larger than n, not linearly (for n < n) 
as in our case. 
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